Sub-Doppler laser cooling of potassium atoms 
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We investigate sub-Doppler laser cooling of bosonic potassium isotopes, whose small hyperfine 
splitting has so far prevented cooling below the Doppler temperature. We find instead that the 
combination of a dark optical molasses scheme that naturally arises in this kind of systems and 
an adiabatic ramping of the laser parameters allows to reach sub-Doppler temperatures for small 
laser detunings. We demonstrate temperatures as low as (25±3)/iK and (47±5)/iK in high-density 
samples of the two isotopes ''^K and ''^K, respectively. Our findings will find application to other 
atomic systems. 
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Sub-Doppler laser cooling of neutral atoms [l| is a key 
technique for the production of ultracold and quantum 
gases. It allows for atoms to be cooled to temperatures 
below the Doppler limit [2[ , not far from the single pho- 
ton recoil energy. This favors the application of further 
cooling techniques, such as Raman or evaporative cool- 
ing, to reach quantum degeneracy. It also realizes a fast 
and effective cooling method for some classes of atomic 
interferometers and clocks [3|. The sub-Doppler cooling 
mechanism arises whenever the atomic ground state has 
an internal structure with state-dependent light shifts. 
Such a situation is typically accompanied by a hyper- 
fine structure of the excited state [J| . The sub-Doppler 
cooling is efficient only if the excited state has a hyperfine 
splitting A either much larger than the natural linewidth 
r, like that for the alkalis Na, Rb and Cs, or smaller than 
r, as for example in Sr [5|. In the intermediate case of 
A '^ r, it can instead be hindered by the presence of 
heating forces or by photon reabsorption [6J . The bosonic 
potassium isotopes fall into this latter category [7|, and 
no efficient sub-Doppler cooling has been observed so far 

We now instead find that sub-Doppler cooling can take 
place also in atoms like K, by employing a near-detuned 
optical molasses and an appropriate strategy to tune the 
cooling laser parameters. We observe that the natural 
depumping towards dark states taking place in this kind 
of systems allows us to reach low temperatures even in 
high-density samples. In experiments on the isotopes 
■^^K and ^^K, we achieve temperatures substantially lower 
than those previously achieved, with an efficiency similar 
to that of most other alkali species. 

In a two-level system, Doppler cooling arises when 
a laser is tuned below the atomic transition frequency, 
where the atoms experience a friction force: F = —av. 
If the presence of Zeeman sub-levels is taken into ac- 
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FIG. 1: Working regions for sub-Doppler cooling of bosonic 
Potassium, a) Level scheme including the relevant hyperfine 
splitting A and the cooling laser detuning 5. b) Calculated 
cooling forces vs the atomic velocity in the various regions 
of (a). Doppler cooling only takes place in regions I and IV, 
while sub-Doppler cooling is active only in regions I and II. 



count [23, a much larger friction arises for small ve- 
locities, leading to temperatures much lower than the 
Doppler limit ksTo = fiY/2. While in principle, the 
lowest achievable sub-Doppler temperatures are indepen- 
dent of the laser detuning 5 [13], the experiments with 
large density samples are performed at large detunings, 
(5 3> r. This requirement arises from the need of keeping 
the scattering rate of photons by individual atoms low, 
in such a way that spontaneously emitted photons may 
not disturb the cooling process [6| . Most atomic systems 
cannot be modeled as simple two-level ones since they 
feature a hyperfine structure like the one in Fig. [T] that is 



relevant for instance to Na, K and Rb. In this case, it is 
commonly thought that 5 must also be smaller than the 
main hyperfine splitting A, since otherwise the presence 
of the other excited states would turn the sub-Doppler 
mechanism into a heating one. As a matter of fact, in 
the case of the bosonic K isotopes, where A Ri 2T — 3T, 
a clear sub-Doppler cooling has not been experimentally 
observed. This can be understood from the nature of the 
optical forces we have calculated for the level structure 
in FigdJ sub-Doppler cooling is active either very close 
to resonance, where heating from photon reabsorption 
might be large, ot ioi S ^ A, where however the veloc- 
ity capture range is very low. Note that in the case of 
r > A, such as for ^^Sr, the sub-Doppler cooling stays 
efficient also for S > A 5]. 

We have now realized that the presence of neighbor- 
ing excited states, however, has also a beneficial effect. 
Indeed it causes a natural depumping of the atomic pop- 
ulation into a dark state, such as the F = 1 ground state 
as shown in Fig. [TJ The atoms can of course be moved 
back into the bright F ~ 2 state by the repumper laser 
but, differently from a pure two-level system, this can be 
done in a controlled way. It is then possible to adjust 
the fraction of atoms in the state coupled to the cooling 
laser in order to optimize the cooling power, while keep- 
ing the reabsorption of spontaneously emitted photons 
under control. Note that the possibility of controlling 
the population of the bright state is absent if A 3> F 
unless an appropriate depumping laser is used [16|. This 
mechanism, which is widely used to trap atomic samples 
at high-density in magneto-optical traps [l7|, turns out 
to be the first essential ingredient for sub-Doppler cooling 
when A « F, since it allows to reach low temperatures 
also when (5 w F. This is apparent from Fig 121 which 
shows the minimum temperature we measured for ^^K 
in near-resonant molasses with a very low intensity of 
the repumping light and a large atomic density (further 
details are given below). The measured temperature is 
well below the Doppler limit already for 5 < F, although 
a further decrease with increasing S is apparent. 

The second important observation is that the sub- 
Doppler cooling survives for detunings larger than the 
Doppler cooling does, as shown for example in the calcu- 
lations of FigUl The lowest temperatures can actually be 
reached only for a range of detunings where the Doppler 
force does no longer provide an efficient cooling. This is 
in principle a problem in experiments, since the veloc- 
ity capture range of the sub-Doppler cooling mechanism 
is usually smaller than the initial thermal velocity, for 
example at the end of the capture stage of a magneto- 
optical trap. Indeed, one normally needs to exploit both 
Doppler and sub-Doppler cooling to achieve low temper- 
atures y. We now find that one can still combine an 
initial Doppler cooling with a final stage of optimal sub- 
Doppler cooling by using a proper dynamical variation of 
detuning and intensity of the cooling laser between the 
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FIG. 2: Optimal sub-Doppler temperatures for '^^K. a) Mea- 
sured temperature without (black squares) and with (red 
dots) the ramping strategy compared to the Doppler theory 
(dashed line). The dotted lines separate the various regions 
as in Fig. H] b) Fraction of atoms remaining in the colder 
component without (black squares) and with (red dots) the 
ramping strategy. 



two regimes of operation. As a matter of fact, the com- 
bination of these two ingredients makes the cooling of K 
as efficient as in the other alkali species. 

We now discuss in detail the experimental strategy. 
The linewidth of the cooling transition for K is F = 
27rx6.0 MHz, which corresponds to a Doppler temper- 
ature To « 145 fjK. The hyperfine splitting A is about 
3.5 F and 2.2 F for ■^^K and "^^K respectively. We perform 
cooling and trapping in a three-dimensional magneto- 
optical trap (MOT) on the D2 transition around 767 nm. 
The trap is loaded with pre-cooled atoms from a two- 
dimensional MOT. After 3 s of loading stage we have ei- 
ther about 2x10"^° atomsof '^^Kor 4x10^ atoms of ''^K at 
temperatures in the 1 mK regime. We then compress the 
cloud via application of a compressed-MOT technique to 
densities around 1 x 10^^ atoms/cm'^. During this initial 
coohng stage we adopt the standard strategy used for 
bosonic potassium [sl, ll3| . We use a detuning larger than 
the whole excited manifold (region IV) with total inten- 
sities as large as 20 Is for both the cooling (Icooi) and 
the repumping (I^ep) beams {I s— 1-75 mW/cm^). This 
allows for a large Doppler capture velocity. Finally, we 
switch off the magnetic field and we cool the cloud in a 
molasses scheme as described below. 

We initially reduce suddenly Irep to 1/100 Icooi and we 
set the repumper beam frequency on resonance with the 
F = 1 — > F' = 2 transition. If we then try to perform 
standard molasses cooling, i.e. by a sudden change of the 
laser parameters to the optimal sub-Doppler cooling val- 
ues, we observe moderate sub-Doppler cooling only for 
small S. A larger S results instead in a bimodal distribu- 
tion of the atomic velocities. A typical instance of such 
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FIG. 3: Sub-Doppler cooling strategy for K. a) Time evo- 
lution of Icooi and S (Irep—l/WOlcooi)- b) Resulting veloc- 
ity distribution measured after a free expansion (black solid) 
compared to that obtained without the ramp (red dashed). 



distribution is shown in Fig. [3}d, as measured by fluores- 
cence imaging after a free expansion of the cloud [ISJ. 
The narrow peak corresponds to sub-Doppler tempera- 
tures, while the broader distribution can be attributed 
to an inefficient Doppler cooling or even to a Doppler 
heating. As shown in Fig. [2] the fraction of atoms in 
the central component decreases as S is increased. A 
more effective strategy consists in first tuning the laser 
to 6 ~ r/2 to provide an initial Doppler cooling and then 
in slowly decreasing the intensity while increasing S, as 
shown in Fig. [3^. This method allows to cool nearly 90% 
of the atoms to lower temperatures, as shown in Fig. [2] 
By minimizing the final temperature we find an optimal 
ramping time of about 10 ms, which corresponds to an 
adiabatic narrowing of the velocity distribution during 
the whole sequence. The minimum temperature attained 
for •^^K is about 25/iK at (5 w 2.5r. It then rises again for 
larger i5, presumably because of the progressive reduction 
of the force. 

We have observed, as shown in Fig. HI that an increase 
of the repumping power prevents the achievement of such 
low temperatures at high density, while the temperature 
does not depend on the repumper power at low density. 
This confirms the role of reabsorption of spontaneously 
emitted photons inside the cloud. 

We have performed analogous measurements on "^^K 
and found a similar behavior. In this case a minimum 
temperature of about 50^K is reached for a detuning 
(5 w F. We have compared the observations with a theo- 
retical estimation of the temperatures achievable for our 
experimental parameters. The optical force, shown in 
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FIG. 4: Measured temperatures for K vs the intensity ra- 
tio of repumping and cooling light, for densities of 4 x 10^° 
atoms/cm'' (black squares) and 8x 10* atoms/cm^ (red dots). 
The heating arising from reabsorption effects at high density 
can be tuned by reducing the repumper intensity. 



Fig.m is calculated from the solution of the optical Bloch 
equations in the semi-classical approximation r, 9|. The 
model is only accurate in ID and for a^ — a^ polariza- 
tions [18|. The agreement with the experimental data 
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FIG. 5: Measured temperatures for '^^K (red triangles) and 
*^K (black squares) and calculated temperatures (lines) vs 
the cooling laser detuning. The dashed line is the predicted 
Doppler temperature. 



shown in Fig. [S] is however rather good. To check the va- 
lidity of these calculations we directly measured the spa- 
tial diffusion coefficient D^ in the optical molasses for ^^K 
[l8| . By a simultaneous measurement of the temperature 
we estimated the friction coefficient as a = ksT/Dx [l|]. 
The magnitude of a « 10^^?ifc^ is in good agreement 
with the calculations. A simple estimation considering 
a two-level system would give a result about 3 orders of 



magnitude larger. We interprete this low friction as a 
result of the macroscopical occupation of the dark state. 
The minimum temperatures at each detuning in Fig. [5] 
were obtained for an almost constant laser intensity. The 
data at small detunings is apparently well described by 
the scaling law observed in several other systems [5|, 120 - 



T = C„ 



nr T I 



Tn 



(1) 



From a combined fit we get C^+^- = 0.20(2) and To — 
9(3)/iK. The C„+„- coefficient is smaller than those mea- 
sured on the two species where sub-Doppler cooling has 
been observed for 5 < F, i.e. *^Er [2^] (0.38(2)) and ^^Sr 
(1.3(3)) but larger than those measured in Rb and 



Cs at large detunings [2l|. The observed scaling for K 



suggests that heating processes due to reabsorption are 
less efficiently suppressed when S is small, as expected. 

To characterize the robustness of the cooling process 
against stray magnetic fields we tried to keep the MOT 
magnetic field on during the cooling procedure. For gra- 
dients larger than 5 G/cm, we reached the Doppler tem- 
perature. This gradient corresponds to an average mag- 
netic field of about 1 G, which is the same characteristic 
value found for the other alkali atoms. 

The techniques described here might be applied to 
other systems, such as the ^S— >^P transitions of ^■^Ca 
H and 1^3 Yb ^, for which A « 3F. Additionally, it 
would be interesting to apply our cooling strategy to Na, 
for which A w 6F. Another interesting case is that of 
an inverted and narrow hyperfine structure like in '^'^K. 
In this case, there are no interfering levels which can 
directly cause heating. However, an increase of S to val- 
ues of the order of A or more leads to a washing out of 
the sub-Doppler cooling mechanism itself, since the de- 
tuning from the various hyperfine levels becomes of the 
same order. We performed numerical simulations of the 
atomic force and we found that this effect leads to a fast 
decrease of the capture velocity with increasing detun- 
ing 18|. This decrease is faster than the one for atoms 



with a large A (or with a large F/A as in the case of 
Sr), and leads to a regime in which the minimum sub- 
Doppler temperature fcsT = Dp/a exceeds the capture 
range. This is presumably the reason for the bimodal 
distribution seen in experiments with "^^K [26[. Weaker 
rates of natural depumping to dark states will possibly 
require forced depumping [16j in high density samples, 
but further experimental and numerical investigations 



are needed [27 1. 



We have shown how the limitations of sub-Doppler 
laser cooling in atomic species with small hyperfine split- 
ting can be overcome by the natural control of the photon 
reabsorption and adiabatic ramping of the laser parame- 
ters. The laser cooling techniques we have developed will 
be easily implemented in all existing experiments with 
potassium atoms. Finally, the direct application of laser 



cooled potassium atoms to interferometric measurements 
might enable a new class of experiments 28 1. 

Note: we recently learned of another report of sub- 
Doppler temperatures in K [29| . 

We gratefully acknowledge contributions by G. Roati, 
S. Ferrari and F. Cataliotti. This work is dedicated to 
the memory of Arturo Bambini who contributed with 
his deep insight to the early investigations on ultracold 
potassium in Florence. This work was supported by 
INFN (MICRA collaboration), by EU (IP AQUTE), by 
ERG (DISQUA and QUPOL grants) and by the ESF and 
GNR (EuroQUASAR program). 



[2] 

[3] 
[4] 



[6] 

[7] 
[8] 

[9] 



[10] 

[11] 
[12] 
[13] 
[14] 
[15] 

[16 

[17] 

[18 
[19 

[20 



Electronic address: llandini@science.unitn.itl 
Present address: Clarendon Laboratory, University of Ox- 
ford, Parks Road, Oxford, OXl 3PU, UK 
S. Chu, Rev. Mod. Phys. 70, 685 (1998); C. Cohen- Tan- 
noudji, ibid. 70, 707 (1998); W. D. Phillips, ibid. 70, 721 
(1998). 

T. W. Hansch and A. L. Schawlow, Opt. Comm. 13, 68 
(1975). 

C. Bord, Metrologia , 39, 435 (2002). 

E. Arimondo, M. Inguscio and P. Violino, Rev. Mod. Phys. 
49, Vol. 1, 31 (1977). 

X. Xu, T. H. Loftus, J. W. Dunn, C. H. Greene, J. L. 
Hall, A. Gallagher, and J. Ye, Phys. Rev. Lett. 90, 193002 
(2003). 

C. J. Cooper, G. Hillenbrand, J. Rink, C. G. Townsend, 
K. Zetie and C. J. Foot, Europhys. Lett. 28, 397 (1994). 

A. Bambini, and A. Agresti, Phys. Rev. A 56, 3040 (1997). 
R. S. Williamson III and T. Walker, J. P, J. Opt. Soc. Am. 
B 121393 1393 (1995). 

C. Fort, A. Bambini, L. Cacciapuoti, F. S. Cataliotti, M. 
Prevedelli, G. M. Tino, and M. Inguscio, Eur. Phys. J. D 
3, 113 (1998). 

M. Prevedelli, F.S. Cataliotti, E.A. Cornell, J.R. Ensher, 
C. Fort, L. Ricci, G.M. Tino, and M. Inguscio, Phys. Rev. 
A 59, 886 (1999). 

L. De Sarlo, P. MaioU, G. Barontini, J. Catani, F. Mi- 
nardi and M. Inguscio, Phys. Rev. A 75, 022715 (2007). 

G. Modugno, G. Ferrari, G. Roati, R. J. Brecha, A. Si- 
moni and M. Inguscio, Science 294, 1320 (2001). 

T. Kishimoto, J. Kobayashi, K. Noda, K. Aikawa, M. 
Ueda, and S. Inouye,Phys. Rev. A 79, 031602(R) (2009). 

J. Dalibard, and C. Cohen- Tannoudji, J. Opt. Soc. Am 

B. 6, 2023 (1989). 
J. Dalibard and Y. Castin, in the Frontiers in Laser Spec- 
troscopy, proceedings of the CXX International School of 
Physics Enrico Fermi edited by T.W. Hansch and M. In- 
guscio (North- HoUand, 1994). 

C. G. Townsend, N. H. Edwards, K. P. Zetie, C. J. 
Cooper, J. Rink, and C. J. Foot, Phys. Rev. A 53, 1702 
(1996). 

W. Ketterle, K. B. Davis, M. A. Joffe, A. Martin and D. 
E. Pritchard, Phys. Rev. Lett., 70, 2253 (1993). 

See the supplementary material for further details. 

T. W. Hodapp, C. Gerz, C. Furtlehner, C. I. Westbrook, 
W. D. Phillips, J. Dalibard, Appl. Phys. B 60, 135 (1995). 

C. Gerz, T.W. Hodapp, P. Jessen, K.M. Jones, W.D. 



Phillips, C.I. Westbrook, and K. Mo.lmer, Europhys. Lett. 

21, 661 (1993). 
[21] Y. Castin, and K. M0lmer, Phys. Rev. Lett. 74, 3772 

(1995). 
[22] A. J. Berglund, S. A. Lee, and J. J. McClelland, Phys. 

Rev. A 76, 053418 (2007). 
[23] R. Maruyama, R. H. Wynar, M. V. Romalis, A. An- 

dalkar, M. D. Swallows, C. E. Pearson, and E. N. Fortson, 

Phys. Rev. A 68, 011403(R) (2003). 
[24] A. Mortensen, J. J. T. Lindballe, I. S. Jensen, P. 

Staanum, D. Voigt, and M. Drewsen, Phis. Rev. A, 69, 

042502 (2004). 
[25] T. Kuwamoto, K. Honda, Y. Takahashi, and T. 

Yabuzaki, Phys. Rev. A 60, R745 (1999). 
[26] G. Modugno, C. Benko, P. Hannaford, G. Roati, and M. 

Inguscio, Phys. Rev. A 60, R3373 (1999). 
[27] Josh W. Dunn and Chris H. Greene, Phys. Rev. A 73, 

033421 (2006). 
[28] G. Varoquaux, R. A. Nyman, R. Geiger, P. Cheinet, A. 

Landragin and P. Bouyer, New. J. Phys. 11 113010 (2009) 
[29] V. Gokhroo, G. Rajalakshmi, R. K. Easwaran and C. S. 

Unnikrishnan, J. Phys. B: At. Mol. Opt. Phys. 44, 115307 

(2011). 



SUPPLEMENTARY MATHERIAL 
TEMPERATURE MEASUREMENT 

Here we give a brief description of our detection strat- 
egy. We image the atomic cloud by florescence imaging 
after time of flight. To do so we release the trap and 
after a variable time we pulse the MOT beams at the 
maximum intensity for 30/is (the high intensity ensures 
the saturation of the atomic transition). We collect the 
scattered light with an imaging system on a CCD cam- 
era. Due to the high optical density of the sample, im- 
ages taken on resonance showed saturation effects and 
hence prevented us from estimating the correct number 
of atoms. Hence we detune the lasers on the blue of 
the F=2— >■ F'=3 transition. Detuning to the red side of 
the transition was avoided due to the proximity of the 
other hyperfine energy levels. To determine the number 
of atoms by this method we first perform spectroscopy on 
the atomic cloud at low density and at the usual inten- 
sity recording the power broadening. To check the valid- 
ity of this technique we also try to perform temperature 
measurements tuning the laser on resonance. Although 
the width of the cloud was affected by saturation effects, 
the resulting temperatures were consistent with the ones 
measured while detuning the laser out of resonance. We 
fit the atomic signal with a 2D Gaussian distribution to 
extract the e~^/^ half- width a of the distribution, mea- 
suring it at different expansion times and henceforth we 
can reconstruct (7{texp), which we fit by the equation 



'^y^exp) 



m:, 



exp' 



(2) 




t (ms) 

exp 



FIG. 6: Typical evolution of the atomic distribution width 
as a function of time after release from the molasses cooling 
beams. The solid line is a fit by Eq. [21 



as well as the initial size of the cloud, which is used in 
order to measure the density of the atomic cloud. 

In order to precisely measure the magnification of our 
imaging system we use two different methods. The first 
one consists of loading the atoms into a magnetic trap 
mounted on a translation stage. The positioning of the 
translation stage has 10/im accuracy. By movement by 
a given amount of the coils and recording the position of 
the center of the cloud on the CCD, we are able to mea- 
sure the magnification, which is found to be: 1/3.23(1). 
The second method consists of allowing the atomic cloud 
to fall due to gravity after switching off the trap and 
recording its position as a function of time. From the fit 
of the position by the law of free fall and assuming for 
the gravitational acceleration: g=9.81m/s^ we find the 
value of the magnification to be 1/3.22(5), which is in 
good agreement with the value determined by the other 
method. 

In case of the bimodal distributions as in Fig. 2 of the 
main paper, we fit just the cold component after a long 
time of flight {texp >10ms). 

DIFFUSION MEASUREMENT 

We measure the spatial diffusion coefficient D^; for 
the Brownian atomic motion during the molasses cool- 
ing phase. This is a useful quantity since it connects 
to the other quantities characterizing the cooling process 
(namely T, a and Dp) via the simple formula: 



D^ 



a 



kuT 



(3) 



From the results of the fit we extract the temperature 



Hence, by measuring both T and Da, one can reconstruct 
the other useful quantities. The measurement is per- 
formed at different final detuning of the molasses cooling 




t(ms) 



FIG. 7: Typical data for the variation of the width of the 
atomic distribution as a function of the time spent inside the 
molasses cooling beams. The solid line is a fit by Eq. O 
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FIG. 8: Measured diffusion coefficients in the optical molasses 
(black squares) compared to the numerically calculated values 
(red dots) as a function of the detuning of the cooling laser. 



beams by allowing the atomic cloud to expand in pres- 
ence of the cooling light for a variable time, and then 
taking an image after 100/zs dark period. By recording 
the variation of the width with time in the molasses cool- 
ing beams, we find it to be in agreement with a diffusion 
process (Fig. [7]). By fitting the evolution of the cloud 
size with the equation: 



ait)' 



oi + 2D,t 



(4) 



we are able to extract the spatial diffusion coefficient. 
From the results in Fig. |8]we see that the agreement with 
the theory is good for our usual experimental parameters 
but the agreement becomes rather poor for smaller de- 
tunings. This might be due to the presence of additional 
heating from rescattered photons since these measure- 



ments were taken at high densities (4x10^° atoms/cm^). 
These values for the diffusion coefficient are a factor 
about of 1000 times higher than the ones measured in 
[l| . The reason for this is the very low population in the 
F=2 level caused by natural depumping and the use of a 
very weak repumping light. 



NUMERICAL SIMULATION 

We perform numerical simulations of the optical Bloch 
equations for multilevel systems to find the optical forces 
acting on the atoms during the molasses cooling. We 
work in ID and in the (t_|_ — U- configuration for the 
laser beams polarizations. This is only an approximation 
to the more complex polarization geometries arising in 
the 3D laser configuration we have in our experiment. 
In practical situations, in a magneto-optical trap, both 
Sisyphus cooling (lin_Llin) and (j+ — a^ polarization gra- 
dient cooling play a role in the cooling process, with the 
former dominating at large detuning and the latter at 
small detuning [3, l3| . In our simulation, only the cooling 
beam is taken into account to find the cooling force. This 
is justified by the very low repumper intensity we use in 
the experiment (about 10~^Is). 

In the determination of the different regions of Fig. 1 of 
the main paper we focus our attention on a velocity range 
corresponding to twice the width of our initial velocity 
distribution (just before the molasses phase the atomic 
temperature is about 1 mK) and we set the cooling laser 
intensity to be Ig. The cooling force was calculated as 
a function of velocity for different detunings. The whole 
range of detunings for which the force was calculated was 
divided into four regions. Region I: This region comprises 
of the detuning values for which the force is opposite to 
velocity for the whole velocity range and hence always 
provides cooling. Region II: The force is opposite to the 
velocity only for small velocities, whereas for higher ve- 
locities the force changes sign, thereby providing heating. 
Region III: In this region, the force provides heating in 
the whole velocity range. Region IV: The force provides 
heating for low velocities, whereas for higher velocities 
the force is opposite to the velocity. In Fig. 1 of the 
main paper we do not report the force behavior for de- 
tuning values in between the F'=2 and F'=0 levels. 

The calculated temperatures of Fig. 5 in the main paper 
are calculated as fc^T = Dp/a where Dp is the momen- 
tum diffusion coefficient and a is the friction coefficient 
calculated as the slope of the force at ti = 0. We estimate 
Dp by using the simple argument for the random step of 
the Brownian motion in momentum space in multilevel 
transitions described in Ref. [4^. 

To compare the efficiency of the cooling process in the 
case of narrow hyperfine structure to the two-level case 
we numerically simulate the cooling force and extract the 
velocity capture range v^, defined as the velocity giving 




capture energy Ec = ^mv"^ for the case of ^^K normalized 
by the 10-fold increased case (black squares), and the in- 
verted case normalized by the same quantity (red dots). 
We show only the value of Ec in the detuning regions in 
which the sub-Doppler cooling mechanism provides cool- 
ing (region I and II of Fig.l in the main paper). This fast 
decrease of the velocity capture range makes the achieve- 
ment of sub-Doppler cooling difficult for optical molasses 
at large detuning for this kind of systems. In conclusion, 
the numerical simulation provides a good understanding 
of the sub-Doppler cooling process and the results ob- 
tained in our experiment. 



FIG. 9: Calculated sub-Doppler capture energy for ^^K (black 
squares) and for a hypothetical atomic species with an in- 
verted hyperfine structure as shown in the inset (red dots). 
Both quantities are normalized to the capture energy of an 
atomic species with A=35r. 



the first local maximum of force at low velocity. We do 
this in three cases: a) a hyperfine level structure like the 
one of '^^K, b) an equal but inverted hyperfine level struc- 
ture, c) the same hyperfine structure but with a 10-fold 
increase in the hyperfine splitting. In FigjS]we plot the 
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